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INTRODUCTION 


Past research with hybrid rockets has suggested that certain motor 
operating conditions are conducive to the formation of pressure oscillations, or 
flow instabilities, within the motor combustion chamber. These combustion- 
related vibrations or pressure oscillations may be encountered in virtually any 
type of rocket motor and typically fall into three frequency ranges: low 
frequency oscillations (0-400 Hz); intermediate frequency oscillations (400-1000 
Hz); and high frequency oscillations (> 1000 Hz) 1 . In general, combustion 
instability is characterized by organized pressure oscillations occurring at 
well-defined intervals with pressure peaks that may maintain themselves, grow, 
or die out. Usually, such peaks exceed ±5% of the mean chamber pressure. 

hybrid motors, these oscillations have been observed to grow to a 
limiting amplitude which may be dependent on factors such as fuel 
characteristics, oxidizer injector characteristics, average chamber pressure, 
oxidizer mass flux, combustion chamber characteristic length L*, and grain 
geometry . Mechanisms which have been postulated as causes of instabilitv in 
hybrid motors include * 

associated with the time lag of vaporization and combustion 
of liquid droplets ' . Inadequate droplet atomization and vaporization, as 
well as inadequate vaporization chamber residence time, are regarded as 
primary contributing factors. These are, in turn, influenced by factors 
such as droplet size and initial injection velocity i.e., injector design), 
and vaporization chamber geometry (e.g., volume, the use of splashb locks , 
etc . ) . This type of instability is characterized by low frequency pressure 
oscillations (5-10 Hz) and it operates generally in the low oxidizer mass 
flow - low chamber pressure regime 4 . 


Instabilities associated with periodic accumulation and break-off of melted 
layers or char layers from the burning propellant surface . Propellant 
formulation and thermal response time are believed to be primary 
contributing factors. This type of instability is characterized by low 
frequency pressure oscillations (1-100 Hz) and occurs with low regression 
rate motor operation . 

Instabilities associated with some boundary layer mechanism which couples 
the combustion response (regression rate) to externally imposed pressure 
oscillations . These instabilities occur in low total mass flux regimes 
(where radiation is the dominant heat transfer mechanism) and in high total 
mass flux and/or low chamber pressure regimes (where finite-rate reaction 
kinetics becomes important) 5 . 

instabilities associated with low impedance LOX feed system design, wherein 
chamber pressure fluctuations result in LOX flow fluctuations. Generally, 
it is believed that fluctuations in LOX flow rate may enhance chamber 
pressure oscillations, but do not initiate them • 


Instabilities associated with shear layer/acoustic interactions in either 
the vaporization chamber or the aft-end mixing chamber of the motor 6 . 

Instabilities associated with the (possibly) time-dependent nature of the 
boundary layer leading edge position (flameholder effects)'. 


PRELIMINARY INSTABILITY AMAT.YSTS 

Low frequency, or so-called ''ngn-acoustic " instabilities have been observed 
in a number of hybrid motor tests*'* 6,2,9,10 # As the name implies, hybrid motors 
share common characteristics with both liquid propellant and solid propellant 
rocket motors. On the other hand, hybrids possess certain unique 
characteristics. For instance, the ability to throttle hybrids over a wide range 
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of operating conditions, coupled with time-varying internal geometry (fuel grain 
configuration) , increases the probability of encountering conditions favorable 
to the onset of instability . In solid propellant motors the combustion process 
is assumed to be localized in a zone close to the propellant surface. The 
combustion is then treated as a boundary condition for the fluid dynamics 
problem, greatly simplifying a stability analysis by isolating chemical reactions 
from the cavity analysis. The combustion dynamics problem is then defined as the 
response of a thin combustion zone at the cavity wall to an imposed flow 
environment. in liquid propellant motors, on the other hand, the flow and 
combustion problems are coupled throughout the entire volume of the combustor. 

Hybrid motor combustion typically occurs in a thin flame sheet located 
within a boundary layer located near the surface of the fuel grain. Fuel is 
partly decomposed and gasified at the solid surface by heat transfer and diffuses 
outward towards the flame zone. At the same time, oxidizer diffuses inward 
towards the flame zone through the boundary layer. At a point where the 
oxidizer-to-fuel concentration (0/F ratio) is slightly on the fuel-rich side of 
stoichiometric, combustion occurs in a zone whose thickness is about 10% of the 
boundary layer thickness. Some hybrid motors also inject oxidizer into a cavity 
located downstream of the solid fuel grain, so that a secondary combustion 
process occurs in this region. Hybrid motors can then be expected to exhibit the 
combustion characteristics of both solid and liquid rocket motors. 

Hybrid Motor L* Instability Analysis 

The approach taken in the present analysis is to develop a modified L* 
instability theory which accounts for the relationship between pressure and 0/F 
ratio in the motor. The analysis proceeds along the lines of that first 
suggested by Beckstead and Price 11 for solid propellant motors, with the 
following simplifying assumptions: 

• The combustion process is assumed to be localized in a zone close to the 
propellant surface. 

• The analysis is for one-dimensional flow of a calorically perfect gas. 
Temperature, molecular weight, and pressure are spacially uniform throughout 
the combustion volume. Molecular weight is either not oscillating or has a 
large oscillation period relative to the flow residence time (quasisteady) • 

• The nozzle discharge is quasisteady, with inlet stagnation conditions 
approximated by static conditions. 

• All flow variables (except the oxidizer mass flow rate) can be expressed as 
the sum of a mean (constant) term and a time-dependent term 

• The oxidizer mass flow is constant. 

The continuity equation applied to the thrust chamber takes the form 


T C d(p/p) + p _ 1 + g _g 

n dt ~P 1+g l+E ci 


( 1 ) 


where p. is pressure, t is time, (-) denotes a mean value, r is the so-called 
motor fill (or flush) time, defined as c 


T 


c 


L*C * 
RT 


( 2 ) 
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and 


a 



(3) 


xs the inverse of the O/F ratio, and n = 1 or y depending on whether the core 
flow xs consxdered to be isothermal or isentropic. In the limit 5 -» oo Eq.(l) 
degenerates to the solid rocket motor relations of Reference 11. Now let 


£ = i + £2 e 9t cos(Pt) (4) 

-§ = x + i e $t CO s(cot+T) (5) 

a a ' ' 


where P is a growth rate constant for oscillations and t is the lead time of the 
fuel-oxxdxzer ratio relative to the pressure. One can define a response function 


♦ - 



( 6 ) 


where a Q and p Q are oscillation amplitudes of the fuel-oxidiaer ratio and 
pressure at some (arbitrary) time t, respectively. Combining Eqs.(l)-(5) and 
solvxng yxelds two conditions (expressed in terms of 0/F ratio) : 


cos (g>t) = 1 


+ 3 ^ 
n 


(7) 




(8) 


Equations (7) and (8) represent a set of 1-D stability oscillator equations for 
hybrid rocket motors. As such, they describe the dynamic behavior of a hybrid 
motor xn terms of the oscillation variables (tt>, P) as functions of the combustion 
dynamics (¥,t), the mean O/F ratio (O/F), and the motor fill time (r >. The 
motor fill time is, in turn, dependent on L . c 

Stabi litv Trends 


If one considers the interval -x < car < vr for the combustion (i.e., fuel~ 
oxxdizer ratio) phase angle, one can quickly show that 


XX— 3 



t > 0 


(9) 


i.e., combustion oscillations lead pressure oscillations for all values of 0/P, 
for either isothermal or isentropic core flow. Considering that the stability 
limit is defined by the condition P = 0, one can show that the oscillator 
equations impose a minimum value constraint on the response function ¥, and that 
this value is a function of the <meanl motor O/F ratio . This constraint is 
simply stated as 


i|r 0 s 1 + 0/F 


( 10 ) 


where the subscript "o" denotes conditions at the motor stability limit. Thus, 
a hybrid propellant motor can be treated (from the standpoint of 1-D if 
stability) as a solid propellant motor with the response function modified to 
account for 0/F ratio. Finally, the expected pressure oscillation frequency can 
be estimated as 


f = 2 £r c ar< = c °S<-^> -f^c.O/F) 


( 11 ) 


so that the expected oscillation frequency depends directly on the system fill 
time t c (or L ) and the average 0/F ratio. 

Sample Results 

Guthrie and Wolf 3 reported non-acoustic instabilities in a series of motors, 
ranging in size from 10,000 lbs thrust to 75,000 lbs thrust, which had been 
developed and tested by the American Rocket Company (AMROC) . They then applied 
solid rocket motor instability theory and compared its predictions with their 
experimental data. In their work, they attempted to correlate data from "several 
solid fuels", and concluded that ¥ « 1.3. According to Eq.(10) o f th e hybrid 
motor theory developed here, the correct value is T = 3.09 (i.e., O/F - 2.09). 
A comparison of predicted vs. observed oscillation frequencies is shown in Figure 
1 . Predicted frequencies are obtained using values of combustor fill time 
reported in Reference 3 . It can be seen that the hybrid motor theory provides 
excellent agreement with the observed values of oscillation frequency. 

SUMMARY AMD CONCLUSIONS 

The hybrid stability theory presented here is a 1-D model derived along the 
same lines as for solid rocket motors, with account taken for a separate oxidizer 
mass flow. It represents a first step towards the understanding of non-acoustic 
oscillations in hybrid rocket motors. Since it is fundamental in nature, it 
necessarily employs a number of simplifying assumptions which restrict the 
applicability and adaptability of the model in more complex situations, and limit 
the type of information which the model can provide. Based upon the limited data 
presented here, the model can predict, to a reasonable degree, the frequency of 
non-acoustic oscillations in simple hybrid motor systems. It must be noted that 
the present theory cannot predict whether or not a given system will oscillate, 
nor can it suggest what measures should be undertaken to reduce or eliminate 
oscillations . 
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Fig. 1: Predicted vs. Observed Oscillation Frequency, AMROC Motors 
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